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1  INTH(  IUTTION 

The  prediction  of  the  effect  of  the  propagation  medium  or.  radio  communi¬ 
cation  systems  operating  in  the  upper  LF  {>100  kHz)  and  MF  (300  kHz  -  3  MHz) 
hands  has  lagged  considerably  behind  predictive  capability  for  the  adjacent 
VLF  (3-30  kHz)  and  Hr  (3  MHz-30  MHz)  bands.  This  has  been  partly  due  to  the 
complexity  of  the  propagation  theory  appropriate  to  the  LF  and  MF  bands  and 
partly  due  to  a  lack  of  operational  requirement.  With  regard  to  the  latter, 
there  has  been  speculation  that  whispering  gallery  type  modes  (i.e.  an  rf  wave 
launched  at  shallow  angles  along  the  lower  boundary  of  the  D-layer  of  the 
ionosphere)  should  be  characterized  by  low  propagation  losses  and  that  the  use 
of  these  modes  could  provide  long  range  air-to-air  communication  links.  '  It 
is  that  speculation  which  lias  motivated  the  present  study.  Although  on  the 
basis  of  the  work  reported  in  references  1  and  2  the  MF  band  appeared  to  be 
the  preferable  of  the  two  bands,  we  have  in  this  study  restricted  the  effort 
to  the  rather  modest  goal  of  waveguide  application  throughout  the  LF  band  for 
a  prototype  daytime  ionosphere. 

Wave  hop'*  techniques  have  been  classically  used  ir.  the  LF  band,  princi¬ 
pally  for  ground-to-ground  transmission.  Some  wave  hop  calculations  for  ele- 

Q 

vated  antennas  have  only  recently  been  reported  and  we  will  make  comparisons 
with  some  of  those  results.  The  present  study  extends  considerably  the  exist¬ 
ing  catalogue  of  daytime-LF  numerical  results.  It  also  demonstrates  the  fea¬ 
sibility  of  using  waveguide  concepts  for  case  studies  of  propagation  under 
daytime  conditions  and,  probably,  for  studies  of  LF  propagation  in  PC  A  and 
nuclear-disturbed  environments.  Use  for  production  applications,  and  possibly 
for  nighttime  conditions,  would  probably  require  improvement  of  the  mr.de 
search  techniques  as  well  as  the  replacement  of  f n 11-  wave  ionospheric 
reflection  calculations  by  either  phase  integral  or  a  n!'in.iMnn  of  phase 


integral  and  full-wave  methods.  One  of  the  attractive  possibilities  m 
future  studies  for  speedier  calculations  in  the  LF  band  and  possibly  ex- 

9-17 

tension  into  the  MF  band  would  be  the  so  called  "hybrid"  methods  “  whereby 
fields  are  calculated  by  combining  waveguide  and  hop  methods. 

In  the  following  section  minor  modifications  of  the  programs  documented 
in  references  13  and  14,  which  allow  for  treatment  of  whispering  gallery 
type  modes,  are  discussed.  In  section  3  the  prototype  ionosphere  is  docu¬ 
mented.  To  check  results  presented  in  reference  8,  a  limited  number  of 
calculations  have  also  been  generated  for  Deeks'^  summertime  profile.  Data 
for  that  profile  are  also  documented  in  section  3.  In  section  4  results  are 
presented  and  following  that  the  conclusions  and  recommendations  arc-  summarized. 


2  SUMMARY  'F 


ON‘S 

The  waveguide  program  documented  in  reference  13  and  the  excitation  ar.i 

height  gain  formulas  documented  in  reference  14  serve  as  the  basis  for  the 

present  study.  Table  1  below  lists  the  real  (  )  and  imaginary  (0.)  part  of 

r  i 

the  eigenangle  and  the  magnitude  of  an  excitation  factor,  |\|,  for  grout  3- 
based  excitation  of  the  first  mode  (mode  numbering  begins  with  the  eigenangle 

having  the  largest  real  part)  as  a  function  of  frequency.  More  complete 

documentation  and  description  of  mode  data  are  presented  in  section  4.  For  the 
present  it  is  only  necessary  to  appreciate  the  fact  that  the  excitation  factor 
gives  a  measure  of  how  well  a  ground-based  antenna  would  excite  the  first 

order  mode  and  that  the  first  order  mode  becomes  more  and  more  earth  detached 

(or  equivalently  more  of  a  whispering  gallery  type  mode)  as  the  frequency 
increases. 


Table  1  Eigenangles  and  Excitation  Factors 


Freq  (kHz) 

9r  (°) 

Si  (°) 

1  X| 

100 

86.098 

-0. 335 

3.2  x  10"9 

150 

86.732 

-0.378 

1  .6  x  10"12 

200 

87.141 

-1.417 

6.3  x  10-17 

250 

87.438 

-- .455 

5.1  x  10-21 

300 

87.670 

-0.493 

2.7  x  10-16 

The  tabulations  in  Table  1  are  consistent  with  the  program  documented  in 
reference  13  and  are  for  the  prototype  daytime  ionosphere  described  in  the 
following  section.  The  significant  feature  of  Table  1  is  that  the  excitation 
factor  does  not  monotonically  decrease  with  frequency  as  it  should  (i.e.  the 
more  earth  -  detached  the  mode  the  poorer  its  excitation  by  a  ground-based 
source).  The  reason  for  the  incorrect  behavior  is  that  the  linear  combination 
of  modified  Hankel  functions  of  order  one  third  used  to  represent  the  height 
gain  at  the  ground  is  in  fact  incorrect  at  the  ground  for  modes  which  are 


highly  earth  detached.  To  correct  for  this  deficiency  we  have  opted  to  throw 


away  earth  effects  altogether  when  the  condition 


Re  |^i  (k/a)(C2  -  aH)3/2/3j  >  12.4 


(1) 


is  met.  The  quantity  within  the  brackets  relates  to  the  degree  of 
evanescence  of  the  height  gain  at  the  ground,  and  the  value  12.4  on  the  right- 
hand  side  of  (1)  has  been  selected  on  the  basis  of  trial  and  error.  It 
requires  the  degree  of  evanescence  at  the  ground  to  be  of  the  order  of  several 
times  10-6  (i.e.  the  magnitude  of  the  modal  height  gain  is  down  by  more  than 
100  dB  from  its  value  near  the  base  of  the  D  layer).  When  condition  (1)  is 
satisfied,  the  plane  wave  reflection  coefficients  uRnd  and  dRdd  referenced  to 
level  d  become 


CHh(qd)  +  F(qd) 
llRHd  CHh(qd)  -  F(qd) 

C  h(q  J  +  i(ct/k)  1/3  h'  (q.) 

R  =  - — -  ,  (3) 

CHh(qd)  -  i(a/k)  /J  h’(qd) 

where  the  subscripts  C  denote  TM  polarization  for  both  the  downgoing  and 
upgoing  wave  and  the  subscripts  I  denote  TE  polarization  for  the  downgoing 
and  upgoing  wave.  Also, 

CR  =  cost  9) 

9  =  eigenangle  referenced  to  level  H 

qz  =  (k/a)2/3  [CH2  +  a(z-H)] 
k  =*  wave  number 
a  =  2/a 

a  =  earth’s  radius 
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d  =  altitude  at  which  modal  equation  is  evaluated 

H  =  altitude  at  which  modified  refractive  index  is  taken  to  be  unity 
h  =  h2  -  exp(i4n/3)h^ 

h2,h1  =  modified  Hankel  functions  of  order  one  third 
i  =  /=J 

F(qd)  =  i  fo.5(a/k)h(qd)  +  (a/k)1/3  h' { qd>j  /n2  ( d) 
n(d)  =  modified  refractive  index  at  height  d 

Rd  =  plane  wave  reflection  coefficient  looking  down  from  level  d 

Tne  subscript  H  on  the  Cs  indicate  that  the  eigenangle  is  referenced  to 
height  H  where  the  modified  refractive  index  is  unity.  Also,  the  prime  on  h 
in  equation  (3)  and  in  the  expression  for  F  denotes  a  derivative  with  respect 
to  the  argument. 

It  should  be  mentioned  that  the  waveguide  program  of  reference  13  does  in 
fact  find  the  correct  eigenangles  for  the  cases  studied  in  this  report  even 

though  the  starting  conditions  at  the  ground  are  incorrect  for  calculation  of 

the  R  values.  The  reason  for  this  is  that  the  admixture  of  incorrect  soiu- 
d 

tions  decays  with  altitude,  and  for  whispering  gallery  type  calculations  the 

mode  equation  is  evaluated  at  the  base  of  the  D  layer  (typically  at  altitudes 

I 

>50  km)  so  that  the  process  is  in  a  sense  self  correcting.  Nevertheless,  it 
seemed  to  us  better  to  use  equations  (2)  and  (3)  under  the  appropriate  condi¬ 
tions  rather  than  risk  the  possibility  of  error  when  parameter  sets  change. 

We  now  turn  to  a  discussion  of  the  mode  sum  calculation  used  in  this 
study.  Apart  from  the  minor  replacements  to  be  discussed,  the  formulas  for 
mode  summing  are  given  in  reference  14.  The  minor  replacements  just  referred 
to  concern  height  gain  replacements  which  must  be  made  when  the  test  given  in 
(1)  is  passed.  Specifically,  equations  <2)  and  (3)  of  reference  14  are  to  he 
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replaced  by 


10 


Finally,  the  height  gain  replacements  (4)  and  (5)  are  also  incorporated 
into  the  waveguide  program  of  reference  13  when  test  1  is  passed.  This  inclu¬ 
sion  eliminates  the  excitation  factor  dilemma  discussed  earlier  in  this  sec¬ 


tion 


3  DESCRIPTION  OF  IONOSPHERIC  PROFILES 

The  majority  of  calculations  were  made  by  using  the  lower  extremities  of  the 

GE-TEMPO  ambient  day  profile.1^  Only  electrons  have  been  included  in  the  pro¬ 
file,  and  their  height  dependence  as  well  as  that  of  the  collision  frequency  are 
shown  in  Table  2.  The  variation  between  tabulated  points  is  assumed  to  be 
exponential.  Though  the  D-region  electron  density  model  has  been  taken  from 
Knapp,  it  has  no  special  significance  and  is  used  in  this  study  simply  to 
demonstrate  the  feasibility  of  carrying  out  waveguide  calculations  throughout 
the  LF  band  for  a  prototype  daytime  ionosphere.  The  profile  has  been  trunca¬ 
ted  at  80  km  since  most  of  the  LF  reflections  occur  below  that  level. 

A  limited  number  of  calculations,  made  to  check  our  results  with  those 
based  on  a  wave  hop  program,  have  been  performed  for  the  Peeks'  summertime  day 
profile  shown  in  Table  3. 

The  GE-TEMPO  and  Peeks'  profiles  are  schematized  in  Figure  1.  Above 
about  64  km  the  GE-TEMPO  number  densities  exceed  those  of  the  Peeks'  profile 
by  factors  varying  between  about  2  and  4.  This  tends  to  make  the  results  for 
the  GE-TEMPO  profile  more  lossy  than  the  Peeks’  profile,  so  that  the  mean  decay 
of  the  signal  for  the  GE-TEMPO  profile  will  be  somewhat  greater.  In  a  compar¬ 
ison  with  data  this  difference  would  be  quite  significant;  but  for  the  princi¬ 
pal  purpose  of  this  study — namely  the  demonstration  of  the  feasibility  of 
carrying  out  waveguide  calculations  throughout  the  LF  band  for  daytime  (and 
presumably  depressed)  conditions — the  difference  is  of  no  significant  conse¬ 
quence. 
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RESULTS 


This  section  contains  a  variety  of  range  and  height  gain  curves  for 
propagation  over  sea.  The  bulk  of  the  curves  are  for  the  GE-TEMPO  profile  and 
are  for  frequencies  of  100,  150,  200,  250  and  300  kHz.  Included  are  rather 

complete  mode  set  tabulations.  The  number  of  modes  range  from  a  dozen  at  100 
kHz  to  28  at  300  kHz.  The  MODESRCH17  algorithm  was  used  to  find  the  mode  set 
for  the  100  and  150  kHz  cases.  Above  150  kHz  however,  numerical  problems  with 
MODESRCH  were  encountered  and  the  modes  were  generated  by  using  the  trace  rou- 
tine  described  by  Ferguson.  Although  the  latter  does  not  infallibly  locate 
all  significant  modes,  missing  modes  can  usually  be  spotted  by  scrutinizing  the 
mode  set  (particularly  as  regards  mode  excitation  and  polarizations  as  well  as 
mode  structure  at  the-  previous  frequency).  They  car.  then  be  located  by  inputs 
of  a  variety  of  prudently  selected  trial  starting  eigenangles .  Obviously  this 
procedure  could  be  made  much  more  cost  effective  by  modifications  to  MODESRCH 
which  would  allow  that  method  to  be  used  at  higher  frequencies.  Very  likely, 
too,  the  trace  routine  of  Ferguson  could  be  improved.  These  are  areas  of 
improvement  recommended  for  future  work. 

Figures  2  through  4  show  comparisons  between  range  calculations  of 

O 

Campbell  and  Jones  ,  who  used  a  wave  hop  analysis,  and  tie  waveguide  calculation 
of  this  study.  Tr.c  cal  ‘illations  are  for  the  Decks'  summertime  day  profile  at 

150  kHz.  Figure  2  applies  to  ground-to-ground  propagation.  Figure  3  is  for 
ground  transmission  to  a  receiver  at  5  km,  and  Figure  4  is  for  ground  trans¬ 
mission  to  a  receiver  at  9  km.  The  waveguide  curves  have  been  generated  for  a 
geomagnetic  field  of  0.4508  Gauss,  a  dip  angle  of  65.82°  and  an  azimuth  of 
0° .  Since  Campbell  and  Jones  specified  only  the  azimuth,  latitude  and 
longitude  for  which  their  calculations  applied,  the  preceding  geomagnetic 
conditions  may  be  at  slight  variance  with  theirs.  Also,  the  digital  input  for 


the  electron  density  and  collision  frequency  may  vary  slightly  from  their 


input.  Differences  of  this  sort  could  be  responsible  for  the  slight  discrep¬ 
ancies  between  the  waveguide  and  wave  hop  calculations.  In  view  of  the  differ¬ 
ence  in  methods,  the  agreement  between  the  two  sets  of  calculations  seems  to 
us  quite  remarkable.  To  our  knowledge  Figures  3  and  4  show  for  the  first  time 
comparisons  between  waveguide  calculations  and  the  wave  hop  method  for  elevated 
antennas. 

Table  4  contains  the  set  of  mode  data  upon  which  the  waveguide  calcula¬ 
tions  in  Figures  2  through  4  are  based.  The  first  column  gives  the  mode  num¬ 
bering  beginning  with  the  eigenangle  which  has  the  largest  real  part  (i.e.  the 
most  grazing  or  earth  detached  mode).  Eighteen  modes  have  been  included  in 
the  set.  The  second  and  third  columns  give  the  real  and  imaginary  parts  of 
the  modal  eigenangle  expressed  in  degrees  and  referenced  to  height  H.  For  the 
calculations  shown  in  Table  4,  H  was  set  to  56  km.  The  fourth  column  gives 
the  modal  attenuation  rates,  which  range  between  about  8.7  dB/Mm  and  36  dB/Mm. 
It  is  interesting  to  note  that  the  modal  attenuation  rate  for  the  ninth  mode 
is  only  about  0.8  dB/Mm  greater  than  the  attenuation  rate  for  the  first  mwde 
(i.e.  the  most  pronounced  whispering  gallery  mode).  There  are  in  this 
instance  about  10  modes  with  attenuation  rates  comparable  to  the  least  atten¬ 
uated  mode  but  no  modes  with  exceptional ly  low  attenuation  rates  as  one  might 
anticipate  for  modes  characterized  by  very  grazing  incidence  angles.  This  is 
because  daytime  ionospheres  are  not  particularly  abrupt  and  there  is  generally 
an  appreciable  ionospheric  absorption  at  altitudes  below  the  height  where  the 
bulk  of  the  reflection  occurs.  We  would  anticipate  this  modal  feature  to  be 
much  the  same  for  PCA  or  artificially  depressed  ionospheres  but  to  be  quite 
different  for  propagation  beneath  an  ambient  nighttime  ionosphere.  The  fifth 
column  gives  the  ratio  of  the  modal  phase  velocity  to  the  speed  of  light  in 
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vuou:~.  The  sixth  and  seventh  columns  give  the  m.i'i:.: :  ..it  ir.d  s  ruse  ( i:i 
f.idia:.  •  of  an  excitation  factor  defined  as 


Me 
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3d 


Vd> 


(7) 


In  the  table,  |A|  is  called  EXTR  MAG  and  |>  is  called  EXTR  ANG.  The  height  was 

taken  to  be  56  km  when  generating  the  results  of  Table  4.  More  precisely,  the 

excitation  factor  in  equation  (7)  is  for  ground-based  vertical  electric  dipole 

excitation  of  the  vertical  electric  field  at  the  ground.  Thus  the  whispering 

gallery  modes  are  weakly  excited  as  are  horizontally  polarized  inodes  relative 

to  vertically  polarized  modes.  The  latter  are  expressed  by  the  polarization 

magnitude  and  angle (in  degrees)  in  the  eighth  and  ninth  columns  of  Table  4.  The 

1  9 

polarizations  are  calculated  by  using  an  equation  given  by  Rapport.  Values 
greater  than  unity  for  the  magnitude  indicate  principally  vertically  polar¬ 
ized  (TM)  waves,  whereas  magnitudes  less  than  unity  indicate  principally 
horizontally  polarized  (TE)  waves.  It  will  be  seen  that  the  first  eight  modes 
contain  comparable  mixtures  of  TE  and  TM  polarization,  "’he  higher  order  modes 
then  divide  into  TE  and  TM  sets.  as  evidenced  by  the  magnitude  of  the 

polarization  as  well  as  the  magnitude  of  the  excitation  factor.  Observe  how 
the  angle  of  the  polarization  alternates  between  a  large  and  small  value.  We 
have  found  this  to  be  a  typical  behavior  throughout  the  LF  hand  and  have  found 
it  useful  in  spotting  potentially  missing  modes.  In  that  regard  the  magnitude 
of  the  excitation  factor  as  well  as  the  magnitude  of  the  polarization  are  also 
quite  useful. 

The  remaining  results  In  this  section  are  for  a  0.6  Gauss  geomagnetic 

field,  an  azimuth  of  45®,  and  a  dip  angle  of  60®. 
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Tables  5  and  7  show  mode  sets 


for  these  geomagnetic  conditions  for  the  Decks'  and  GE-TEMPO  profiles.  The 
tables  apply  to  150  kHz  and,  in  particular.  Table  5  applies  to  the  Deeks' 
summertime  day  profile  and  Table  7  to  the  GE-TEMPO  day  profile.  In  each 
instance  eighteen  modes  are  shown.  For  the  Deeks'  profile  the  attenuation 
rates  range  between  about  8.7  dB/Mm  and  36  dB/Mm  (just  as  they  did  for  the 
geomagnetic  conditions  which  apply  to  Table  4)  whereas  they  range  between 
about  11  dB/Mm  and  46  dB/Mm  for  the  GE-TEMPO  profile.  The  additional  loss  for 
the  GE-TEMPO  profile,  as  explained  previously,  would  be  anticipated  on  the 
basis  of  its  greater  ionization  above  64  km  (see  Figure  1).  The  consequences 
of  the  higher  attenuation  rates  on  mode  sum  plots  are  shown  in  Figures  5 
through  7.  Those  curves  show  range  plot  comparisons  for  the  Deeks'  and  GE- 
TEMPO  profiles.  Figure  5  is  for  a  ground  based  transmitter  and  receiver. 
Figure  6  is  for  transmitter  and  receiver  altitudes  of  30  km  and  Figure  7  is 
for  transmitter-receiver  altitudes  of  50  km.  At  the  more  distant  ranges 
differences  in  excess  of  10  dB  occur.  The  thing  most  amazing  to  us,  and 
unexplained,  is  the  coincidence  in  the  null  and  maxima  locations  for  the  two 
rather  disparate  profiles.  The  proliferation  of  mode  structure  with  increasing 
altitude  of  the  transmitter  and  receiver  is  evident,  though  the  deepest  nulls 
rather  surprisingly  occur  for  the  transmitter-receiver  altitude  of  30  km. 
Picking  up  more  of  the  whispering  gallery  modes  with  the  transmitter-receiver 
combination  at  50  km  tends  to  fill  in  the  deep  nulls  occurring  with  the  30  km 
combination. 

All  of  the  remaining  results  in  this  section  are  for  the  GE-TEMPO  day 
profile.  Tables  6  through  10  give  the  mode  data  sets  at  100,  150,  200,  250 

and  300  kHz.  The  number  of  modes  range  from  a  doz-'-1  at  100  kHz  to  28  at  300 
kHz.  Minimum  modal  attenuation  rates  range  from  about  7.7  dB/Mm  at  100  kHz  to 
about  19  dB/Mm  at  300  kHz.  Mode  spacing  for  the  real  part  of  the  eigenangle 


is  floral ly  comparable  to  or  greater  than  a  fex  hundredths  of  a  degree  and 
for  the  imaginary  part  of  the  eigenangle  the  mode  spacing  is  generally  compar¬ 
able  to  or  greater  than  a  few  thousandths  of  a  degree.  There  can  be  excep¬ 
tions  however.  The  seventeenth  and  eighteenth  modes  for  the  300  kHz  case 
(Table  10)  have  real  parts  which  differ  by  only  0.004°  and  have  imaginary 
parts  which  are  identical  to  the  number  of  places  printed  out  (i.e.  to  a 
thousandth  of  a  degree) .  This  points  out  the  rigid  demands  on  the  program 
from  the  standpoint  of  eigenvalue  resolution.  It  will  also  be  noted  from  the 
tabulations  that  the  excitation  factors  for  the  whispering  gallery  type  inodes 
show  a  monotonic  decrease,  unlike  the  behavior  in  Table  1,  with  frequency. 

Figures  8  through  12  show  range  plots  for  frequencies  of  100,  150,  200, 
250  and  300  kHz.  On  each  plot  are  curves  for  ground-to-ground  transmission, 
for  transmitter  and  receiver  at  30  km  and  for  transmitter  and  receiver  at  50 
km.  The  null  in  the  ground-to-ground  transmission  curves  which,  depending 
upon  frequency,  falls  in  the  range  from  about  800  to  1400  km  is  a  manifesta¬ 
tion  of  the  ground  wave  and  first  hop  sky  wave  interference  null.  Even  up 
through  300  kHz  the  ground-to-ground  curves  show  relatively  little  modal 
structure, indicating  that  only  a  few  modes  are  require)  for  that  configura¬ 
tion.  The  mode  structure  is  considerably  more  complicated  for  the  elevated 
transmitter  and  receiver  cases.  On  the  basis  of  the  mode  picture  this  would 
be  anticipated  since  whispering  gallery  type  modes  tend  to  become  more  influ¬ 
ential  with  terminal  elevation.  On  the  basis  of  a  wave  hop  picture  the  added 
s*  r  . -turn  would  be  anticipated  since  more  multi)  lath  jx/suil.-i  1  it  le-.  ex  i 
when  the  terminals  are  elevated.  When  the  *-ransmitter  and  receiver  are 
or.  the  ground,  for  exam;  1<-,  only  one  path  applies  to  a  smile  i onospher  1  r 
reflection.  When  the  terminals  are  elevated,  on  t  he  *  '  '.and,  there  exist 
four  paths,  or  hops,  linking  transmitter  and  receiver  which  correspond  to  a 


single  ionospheric  bounce.  Though  the  absolute  signal  levels  and  the  location 


of  the  nulls  would  be  quite  sensitive  to  the  ionospheric  model,  the  severity 
of  multimode  interference,  depth  of  nulls,  etc  is  probably  quite  realistical¬ 
ly  modeled  by  the  results  shown  in  Figures  8  through  12.  Generally  it  will  be 
observed  that  the  deepest  fades  occur  for  the  transmitter  and  receiver 
altitudes  of  30  km.  For  ground-to-ground  transmission  the  fact  that 
whispering  gallery  inodes  do  not  play  a  role  reduces  the  modal  interference  but 
at  the  same  time  also  reduces  the  mean  signal  level.  The  50  km  to  50  km 
transmission  link  where  the  whispering  gallery  modes  play  their  fullest  role 
shows  the  largest  signal  strengths  on  the  mean;  but  as  the  figures  show, 
rather  deep  fades  can  be  expected  for  that  configuration  as  well. 

Figures  13  through  17  show  the  height  behavior  of  the  total  field  for 
frequencies  of  100  kHz,  150  kHz,  200  kHz,  250  kHz  and  300  kHz.  Each  figure 
contains  three  curves.  One  is  for  a  ground  based  transmitter,  another  is  for 
a  transmitter  at  30  km  and  the  third  is  for  a  transmitter  at  50  km.  All 
curves  are  for  a  range  of  2  Mm.  The  first  thing  that  is  striking  about  these 
plots  is  the  depth  of  the  nulls  that  would  be  expected  for  the  elevated  trans¬ 
mitter  cases.  The  altitude  location  of  these  nulls  would  be  sensitive  to  the 
range  as  would  be  the  absolute  field  strength  at  any  given  altitude.  For  the 
case  selected  (i.e.  a  range  of  2  Mm),  it  is  true  that  for  all  frequencies  the 
strongest  ground  signal  obtains  with  the  ground  based  transmitter  and  that  the 
strongest  signal  at  50  km  results  when  the  transmitter  is  at  50  km.  It  is 
true  also,  for  the  cases  examined,  that  the  ground  based  transmitter  yields 
the  lowest  3ignal  level  at  50  km. 

In  addition  to  speedier  techniques  for  determining  the  ionospheric 
reflection  coefficients  as  well  as  more  automated  mode  search  capability,  LF 
propagation  may  be  better  treated,  especially  under  nighttime  conditions,  by 
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alternative  methods.  The  wave  hop  method  already  mentioned 


one  such  possi¬ 


bility.  Another  possibility  would  be  something  akin  to  the  hybrid  method  of 
Felsen  and  coworkers,  where  a  mixture  of  modes  and  hops  are  used  to  describe 
the  field.  The  remainder  of  the  section  has  been  structure!  with  that  possi¬ 
bility  in  mind.  In  particular  we  will  suggest  a  possible  scheme  for  deciding 
on  the  modes  to  be  included  though  we  will  leave  for  future  study  the  treat¬ 
ment  and  inclusion  of  hops.  For  simplicity  we  also  restrict  the  discussion  to 
range  considerations  with  the  transmitter  at  the  same  altitude  as  the 
receiver. 

Modes  which  are  highly  evanescent  at  the  terminal  locations  contribute 
little  to  the  mode  sums.  The  degree  of  evanescence  is  determined  roughly  by 
the  factor  exp|-2/3( k/a) (a  where  Az  is  measured  upwards  from  the 
terminal  altitude.  For  the  sake  of  argument  let  us  say  we  require  the  degree 
of  evanescence  to  be  20  dB.  Then  we  require  roughly  that 

8.68  [|  (k/a)  (a  Az) 3/2 ]  =  20  .  (7) 
Taking  300  kHz  for  sample  calculations,  equation  (7)  gives 


Az  -  9.88  km  . 


(8) 


The  condition  that  the  wave  just  evolves  into  an  evanescent  stage  at  zn  +  A 
where  zQ  is  the  transmitter-receiver  altitude,  is 


2 

C, ,  +  a  ( z  +  Az  -  H)  =  0  . 
H  o 


(9) 
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The  imaginary  part  of  is  ignored  in  these  estimates. 


From  equation  (9)  we 


find  for  zQ  *  0 ,  30,  and  50  km  the  following  (H  was  taken  to  be  55  km  for  all 
calculations  involving  the  GE-TEMPO  day  profile): 

zQ  =  0  9r  =  83.16° 

zo  =30  0r  =  86.05°  (10) 

zD  =  50  9r  =  90°  . 

The  last  of  equation  (10)  is  the  interpretation  for  the  real  part  of  9  when 
2 

Cjj  given  by  equation  (9)  is  negative.  Equations  (10)  give  the  upper  bound  on 

the  eigenangle  search  for  the  three  transmitter  receiver  configurations.  It 

is  suggested  that  the  lower  bound  be  selected  as  follows:  The  hop  calcula- 

\ 

tions  are  most  easily  implemented  when  asymptotic  expansions  can  be  used  for 
the  modified  Hankel  functions  of  order  one  third.  We  thus  determine  the  lower 
bound  on  CH  by  requiring  that 


<k/a)2/3  Jc2  +  a  ( -  H)]  = 


(11) 


The  left-hand  side  of  equation  (11)  is  the  argument  of  the  Hankel  functions  of 
order  one  third,  and  the  number  5  on  the  right  hand  side  has,  for  the  sake  of 
example,  been  assumed  sufficiently  large  relative  to  one  to  justify  asymptotic 
expansions  of  the  Hankel  functions  of  order  one  third.  Again  the  imaginary 
part  of  Cjj  is  ignored  in  these  estimates.  From  equation  (11)  we  find 
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0 


31  .03 


zQ  =30  Sr  =  83.05°  (  12 

z0  =50  3r  =  84.76°  . 

Combining  the  results  from  (10)  and  >.±2)  ,  the  ranges  over  which  modes  would  be 
located  are 

zQ  =  0:  91.08°  <  9r  <  83.16°:  Modes  17-28 
zQ  =  30:  83.05°  *  5r  <  86.05°:  Modes  5-18 
zQ  =  50:  84.76°  <  3r  <  90°;  Modes  1-3  - 

Calculations  based  on  these  combinations  are  shown  in  Figures  13  through  20, 
where  they  are  compared  with  the  results  of  the  full  nwde  set. 

The  results  for  the  limited  mode  set  for  the  zQ  =  0  case  are  indistin¬ 
guishable  from  the  results  for  the  full  modes  set.  This  is  because  the  lower 
limit  81.08“  on  9r  is  within  the  region  of  very  lossy  modes.  As  a  matter  of 
fact  even  the  restricted  mode  range  17-28  includes  many  rtore  n»des  than  are 
necessary  to  adequately  represent  the  zQ  =  0  case.  The  disparity  between  the 
full  and  abbreviated  mode  set  results  for  the  cases  zQ  =  30  and  50  km  points 
out  the  importance  of  higher  order  modes  or  hop  contributions.  Supplementing 
the  restricted  mode  calculations  with  hops  via  the  methods  of  Felsen  and  co¬ 
workers  is  an  interesting  possibility. 
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CONCLUSIONS 


Waveguide  calculations  have  been  carried  out  throughout  zhe  LF  band  for 
a  daytime  ionospheric  model.  On  the  basis  of  the  study  it  would  appear  feasible 
at  the  present  time  to  perform  case  studies  for  additional  daytime  models  as 
well  as  for  PCA  or  artif ically  depressed  ionospheres.  Production  runs,  how¬ 
ever,  would  depend  upon  the  outcome  of  two  developments:  !i)  Replacement  of 
full-wave  ionospheric  reflection  calculations  by  either  phase  integral  methods 

or  by  a  combination  of  phase  integral  and  full-wave  methods.  (2)  Improve- 

1  *7 

merit  and  automation  of  mode  search  techniques.  MODES RCH  was  found  to  have 

numerical  difficulties  at  frequencies  above  about  ISO  kHz,  and  the  TRACE  rou- 
18 

time  of  Ferguson  can  switch  without  warning  from  tracing  one  factor  of  the 

modal  equation  to  tracing  the  other.  Very  likely,  too,  the  TRACE  routine 

would  totally  miss  modes  over  low  conductivity  terrain.  Extension  of  cal cu la t iona 1 

capability  for  nighttime  ionospheres  and  into  the  MF  band  for  both  day  and 

3  -7 

night  ionospheres  might  require  alternative  methods.  The  wave  hop  and 
9-12 

"hybrid"  methods  are  two  such  possibilities. 

Mode  sums  and  height  gains  generated  in  this  study  point  out  the  severity 
of  fading  to  be  expected,  particularly  in  the  upper  LF  band,  when  elevated 
transmitters  and  receivers  are  involved.  For  the  prototype  daytime  ionos¬ 
pheres  considered  in  this  study,  the  number  of  modes  ranged  from  a  dozen  at 
100  kHz  to  28  at  300  kHz.  The  results  generated  here  should  form  a  useful 
basis  against  which  to  compare  results  of  alternative  methods. 

RECOMMENDATIONS 

1.  Improve  mode  search  methods  and  simplify  ionospheric  full  wave  re¬ 
flection  calculations  for  application  of  wavequide  methods  ?•  the  upper  LI' 
band  and  above . 

2.  Investigate  the  utility  of  "hybrid"  methods  wh>-r* by  f :>  ids  are  cal¬ 
culated  by  combining  waveguide  and  wave  hop.  methods. 
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GE-TEMPO  Daytime  Electron  Density  (Ne)  and 
Collision  Frequency  (v  )  Profiles 


.tude  (km) 

Ne(cm"3) 

v0(s  1) 

SO 

4.2'’  x  10" 1 

1.29  x  108 

55 

6 . 50 

60 

8.43  x  101 

65 

2.07  x  102 

70 

5 . 72  x  102 

75 

1  .34  x  103 

80 

2.30  x  103 

9.99  x  105 

23 


Deeks'  Summertime  Day  Electron  Density  (Ne)  and 
Collision  Frequency  (vg)  Profiles 


Altitude 

(km) 

90  .00 

89.50 
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88.00 
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86.00 
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Mg ( cm” 3 ) 
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4.60+002 
4.40+002 
4.15+002 
3  .90+002 
3.60+002 
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ve(s  1  ) 

2  x  105 
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(km) 
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66.50 

66.00 

65.50 
65.00 

64.50 
64.00 

63.50 
63.00 
62  .50 
62.00 
61  .50 
61  .00 
60  .50 
60.00 

59.50 
59.00 

58.50 
58.00 
57  .50 
57.00 

56.50 
56.00 


Ne ( cm" 3  * 


2.90+002 
2.45+002 
1  .90+002 
1 .60+002 
1  .35+002 
1  .20+002 
1  .10+002 
1  .00+002 
9.80+001 
9.60+001 
9.70+001 
1  .00+002 
1  .10+002 
1  .20+002 
1  .35+002 
1  .45+002 
1  .60+002 
1 .65+002 
1  .75+002 
1  .75+002 
1  .70+002 
1 .65+002 
1  .55+002 
1  .40+002 
1  .20+002 
1.05+002 
8.30+001 
6.40+001 
4  .60+001 
3.00+001 
1  .85+001 
1  .00+001 
5.44+000 
2.96+000 
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Figure  1.  GE-TEMPO  (dashed)  and  Deeks'  summertime  (solid)  electron 
density  profiles. 


Azunuth  =  0° ,  dip  -  65. 82°,  B  =  0.4508  gauss 
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Figure  7.  Waveguide  comparisons  for  Deeks '  (dashed)  and  GE-TEMPO  (solid) 
profiles  (ISO  kHz).  Transmitter  and  receiver  altitudes  =  50  km. 
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Figure  8.  Range  plots  for  GE-TEMPO  day  (100  kHz). 


G.E. -TEMPO  DRY 
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G.E. -TEMPO  DRY 
FREQ  -  200  KHZ 
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Figure  10.  Range  plots  for  GE-TEMPO  day  f200  kHz). 


Figure  11.  Range  plots  for  GE-TEMPO  day  (250  kHz). 
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Figure  14.  E  as  a  function  of  height  for  GE -TEMPO  day  (150  kHz) . 
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Figure  15.  E  as  a  function  of  height  for  GE-TEMPO  day  (200  kHz) . 
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Figure  16.  E  as  a  function  of  height  for  GE-TEMPO  day  (250  kHz). 
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FIELD  IdB  ABOVE  1  LlV/m  FOR  1  kW) 


GE  -  TEMPO  DAY 
FREQ  =  300  kHz 


Figure  18.  Range  plot  comparisons  between  full  and  truncated  mode 
sums  (300  kHz).  Transmitter  and  receiver  altitudes  =  0  km. 


Figure  19.  Range  plot  comparisons  between  full  and  truncated  mode 
sums  (300  kHz).  Transmitter  and  receiver  altitudes  =  30  km. 
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Figure  20.  Range  plot  comparisons  between  full  and  truncated  mode 
sums  (300  kHz) .  Transmitter  and  receiver  altitudes  =  50  km. 
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